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Abstract 



This paper studies one-dimensional Ornstein-Uhlenbeck processes, with the distinguish- 
ing feature that they are reflected on a single boundary (put at level 0) or two boundaries 
(put at levels and d > 0). In the literature they are referred to as reflected OU (ROU) and 
P^ ' doubly-reflected OU (DROU) respectively. For both cases, we explicitly determine the decay 

Ph . rates of the (transient) probability to reach a given extreme level. The methodology relies 

r^ ' on sample-path large deviations, so that we also identify the associated most likely paths, 

j^ ■ For DROU, we also consider the 'idleness process' Lt and the 'loss process' Ut, which are 

the minimal nondecreasing processes which make the OU process remain ^ and ^ d, 
respectively. We derive central limit theorems (CLTs) for Ut and Lt, using techniques from 
stochastic integration and the martingale CLT. 
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1 Introduction 

Ornstein-Uhlenbeck (OU) processes are Markovian, mean reverting Gaussian processes. They 
well describe various real-life phenomena, and allow a relatively high degree of analytical 
tractability As a result, they have found wide-spread use in a broad range of application do- 
mains, such as finance, life sciences, and operations research. In many situations, though, the 
stochastic process involved is not allow to cross a certain boundary, or is even supposed to 
remain within two boundaries. The resulting reflected (denoted in the sequel by ROU) and 
doubly-reflected (drou) OU processes have hardly been studied, though, a notable exception 
being the works by Ward and Glynn IT3l[14l [T5l, where ROU processes are used to approximate 
the number-in-system processes in M/M/1 and GI/GI/1 queues with reneging under a spe- 
cific, reasonable scaling; the DROU process can be seen as an approximation of the associated 
finite-buffer queue. Srikant and Whitt ||TT| also show^ that the number-in-system process in a 
GI/M/n loss model can be approximated by ROU. For other applications, we refer to e.g. the 
introduction of [7] and references therein. 

As know^n, the OU process is defined as the unique strong solution to the stochastic differential 
equation (SDE): 

dXt = {a- 'yXt)dt + adBt, Xq = x E R, 

where a S R, 7, cr > and Bt is a standard Brownian motion. This process is mean-reverting 
towards the value a/7. To incorporate reflection at a lower boundary 0, thus constructing ROU, 
the following SDE is used, where we throughout the paper additionally assume a > 0, 

dYt = (a - -iYt)dt + adBt + dLt, Yo = x ;^ 0, 

where Lt could be interpreted as an 'idleness process'. More precisely, Lt is defined as the 
minimal nondecreasing process such that Yt ^ for t ^ 0; it holds that L ^i l^Yt>o}'^Lt = 
for any T > 0. 
Likewise, reflection at two boundaries can be constructed. DROU is defined through the SDE 

dZt = {a- -fZt)dt + adBt + dLt - dUt, Zq = x £ [0, d], 

where Ut is the 'loss process' at the boundary d, i.e., we have Tq y, l^Zt>o}'^Lt = as well as 
L ^1 l^Zt<d}'^Ut = for any T > 0. In the case of DROU we assume that the upper boundary 
d is larger than a/7 throughout this paper, to guarantee that hitting d does not happen too 
frequently (which is a reasonable assumption for most of applications). For the existence of 
a unique solution to the above SDEs with reflecting boundaries, w^e refer to e.g. lfT2| . In the 
context of queues with finite-capacity, Ut is the continuous analog to the cumulative amount 
of loss over [0, t], and that explains why we refer to it as the 'loss process'. 

A first objective of this paper is to obtain insight into transient rare-event probabilities. We do 
so for an ROU process with 'small perturbations', that is, a process given through the SDE 

dy/ = (a - 7y/)di + V^adBt + dL^ 



with e > typically small. The transient distribution (at time T ^ 0, for any initial value 
x ^ 0) of the OU process being explicitly known (it actually has a Normal distribution), we 
lack such results for the ROU process. (As an aside, we note that the stationary distribution of 
ROU is known |14|; it is a truncated Normal distribution.) The above motivates the interest in 
large-deviations asymptotics of the type 

lim e log P(y^ ;?b\Y^ = x), (1) 

for X ^ 0, r ^ 0, and b > E{Y^ \ Y^ = x) (so that the event under consideration is rare). We 
follow the method used for computing blocking probabilities of the Erlang queue in [lOJ, that 
is, relying on sample-path large deviations. In our strategy, a first step is to study the above 
decay rate for the 'normal' (that is, non-reflected) OU process. This decay rate is computed 
as the solution of a certain variational problem, relying on standard calculus-of-variations: 
it minimizes an 'action functional' over all paths / such that /(O) = x and f{T) ^ b. The 
optimizing path /* has the informal interpretation of 'most likely path' (or: 'minimal cost 
path'): given the rare event under study happens, w^ith overwhelming probability it does so 
through a path 'close to' /*. The next step is to observe that /* does not hit level between 
and T, and hence the decay rate that we found also applies for ROU (rather than OU). 
The computations for OU are presented in Section |2l The results are in line with what could 
be computed from the explicitly known distribution of X^ conditional on Xq = x, but pro- 
vide us in addition w^ith the most likely path. Section |3] then focuses on the computation of 
the decay rate for ROU. Above we described the intuitively appealing approach we followed, 
but it should be emphasized that at the technical level there are some non-trivial steps to be 
taken. The primary complication is that the local large-deviations rate function at the reflecting 
boundary is different from this function in the interior [3J. Inspired by (H, w^e derive explicit 
expressions of the large-deviations rate function for ROU by properties of the reflection map in 
the deterministic Skorokhod problem. Unfortunately, calculus-of-variation techniques cannot 
be used immediately to identify the most likely path; this is due to the fact that we need to min- 
imize over all non-negative continuous paths. However, the non-negativity of the optimizing 
path for the OU process facilitates the computation of the decay rates for ROU. In Section 4, we 
compute the decay rate for DROU by the same strategy as the one for ROU. 

The second part of the paper focuses on DROU, with emphasis on properties of the loss process 
Ut (and also the idleness process Lt), for t large. Zhang and Glynn's martingale approach, as 
developed in |16|, is employed to tackle a problem of this type. With h{-) being a twice contin- 
uously differentiable real function, we apply Ito's formula on h{Zt) and require h{-) to satisfy 
certain ordinary differential equations (ODE s) and specific initial and boundary conditions in 
order to construct martingales related to Ut and Lt. The presence of Zt in the drift term leads 
to ODE s with nonconstant coefficients, which seriously complicates the derivation of exact so- 
lutions. In Section|5]we use this approach to identify a CLT for Ut- we find explicit expressions 
for qu and i]ij such that {Ut — qut)/Vt converges to a Normal random variable with mean and 
variance ry^; a similar result is established for Lt. In Section [6l we discuss the corresponding 



large deviations probability, that is the probability that Ut/t exceeds a given threshold c larger 
than Qu, for t large. 

2 Transient asymptotics for Ornstein-Uhlenbeck 

The primary goal of this section is to compute the decay rate Q with V^ replaced by X"^; 
in other words, we now consider the OU case (that is, no reflection). Before we attack this 
problem, we first identify the OU process' average behavior. To this end, we first describe the 
so-called 'zeroth-order approximation' of one-dimensional diffusion processes. The SDE (more 
general than the one defining OU) we here consider is 

dJl = b{Jl)dt + yfe(T{Jl)dBt, Jq=x, 

and the corresponding ODE is 

dx(t) = 6(x(t))dt, x(0) = X. 

Theorem 2.1 IS Thm. 2.1.2] Suppose that b{-) and a{-) are Lipschitz continuous and increase no 
faster than linearly, i.e., 

[b{x) - h{y)]^ + [a{x) - a{y)f ^ K^\x - y\\ 

b^(x) + a'^{x) ^/C2(l + |a;|2), 
where K is a constant. Then for allt > and e > Owe have 

^Jl -x{t)\^ ^ ea(t), 

where a{t) is a monotone increasing function, which is expressed in terms of\x\ and K. Moreover, for 
allt > and 5 > 



limP sup I J| - x{s)\ > 6] =0. 
<:-^0 \Os:ss£t / 

In the specific case of OU processes, the corresponding small perturbation process X| (on a 
finite time interval) satisfies 

dX^ = {a- -fX^)dt + VeadBt, Xq = x ^ 0. (2) 

It is readily checked that the limiting process x{t) is given by 

x(t) = a — 7x(t), x{0) = X, 

which has the solution 



Note that x{t) = KX^. Popularly, as e | 0, with high probability X^ is contained in any 6- 
neighborhood of x{t) on the interval [0, T]. Assuming that b > x{T), it is now seen that the 
probability of our interest, of which we wish to identify the decay rate, relates to a rare event. 

We now recall the Freidlin-Wentzell theorem |2, Thm. 5.6.7], which is the cornerstone behind 
the results of this section. To this end, we first define Cp^rj (M) as the space of continuous 
functions from [0,T] to M, with the uniform norm ||/||oo := supjgjo,T] 1/(01 ^^id the metric 
d{f, g) '■= 11/ — ^lloo- The Freidlin-Wentzell result now states that X^ satisfies the sample-path 
large deviations principle (ldp) with the good rate function 

j^(y) .= I (2^^')"' /o'(/'W - « + If it)? At if / € H,, 

where H^ := {/ : f{t) = a; + /g 0(s)ds, cp G L2([0, T])}. The LDP states that for any closed set F 

and open set G in (C[o,r](J^)) II ■ I loo)/ 

limsupelogP(Xi' G F) ^ - inf I^{f), 
lim inf e log F{X^ £ G) ^ - inf 4 (/) . 

These upper and lower bounds obviously match for /^^ -continuity sets S, that is, sets S such 

thatinf/gci5 4(/) = inf/gints'4(/)- 

We now return to the decay rate under consideration. Let us first introduce some notation, 
foUow^ing standard conventions in Markov process theory. We w^rite P^; {E) for the probability 
of an event E in terms of the process X'^ if this process starts in x. We will mainly w^ork with 
a fixed time horizon T > and write X_ for {Xt, t G [0, T]}. Our first step is to express the 
probability under study in terms of probabilities featuring in the sample-path LDP. Observe 
that we can write F^{X^ ^ b) = P^(Xf G 5), with 

S:=\J Sa, Sa := {/ G qo,T] W : /(O) = x, f{T) = a} . 

Below we first solve a cakulus-of-variation problem to find inij^Sa ^x{f) explicitly. Secondly, 
we prove that S is an /a^. -continuity set. A combination of these findings gives us an expression 
for the decay rate. 
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Moreover, f*{t) ^ Oont £ [0, oo) when the starting point x ^ 0; f*{t) G [0, d]ont G [0, T] when the 
starting point x G [0, d], a G [0, d] and 0/7 < d. 



Proof Obviously, 

m|^ Uf) = inf {^l U'it) + lf{t) - afdt, / G //. n S,| . 

According to Euler's necessary condition flOl Thm. C.13], the initial condition and the bound- 
ary condition, we have that the optimizing path satisfies 

/"(0-7'/(i) + «7 = 0, /(0) = x, f{T) = a. 

The general solution of the ODE (unique up to the choice of the two constants) reads 

/(t) = C7ie^* + C72e-^* + -. 

7 

It is now readily checked that the stated expression follows, by imposing the initial condition 
and the boundary condition. Hence, 

We proceed with proving that /*(t) ^ and /*(i) G [0, d] on t € [0, oo) under the two stipu- 
lated assumptions. First we note that x{t) = xe""^* + (1 — e~'>'*) a/7, a convex combination of 
X and a/7. Since both of these are nonnegative by assumption, so is x{t). For f*{t) we have 
the following alternative expressions with q{t) := sinh(7t)/sinh(7T), as a direct computation 
shows: 

nt)=x{t) + {a-x{T))q{t) 

= q{t)a + (e-^* - g(t)e-^^)x + (l - e"^* - q{t){l - e"^^)) -. 

7 

It follows from the first equality that /*(t) ^ x{t), because a ^ x{T), and hence f*{t) is non- 
negative. Moreover, the second equality shows that f*{t) is a convex combination of a, x and 
a/7, see below. Since all three of these are assumed to be less than d, the same holds true for 
f*{t). Finally we show that we indeed have the claimed convex combination, by showing that 
all coefficients are nonnegative and sum to one. The latter is obvious, as well as q{t) G [0, 1]. 
Furthermore e~^* — q{t)e~^'^ ^ (1 — q{t))e~'^'^ ^ 0. To prove that the third coefficient is non- 
negative we use the basic equality 

. , , , (l + e")(l-e-") 
smh(x) = ^ — -. 



Then observe that 

This completes the proof. D 



Proposition 2.3 S is an Ix-continuity set. 

Proof Consider the topological space (C[o_t](I^))''")/ where the topology r is induced by the 
metric d{f,g). We next consider Sx = {f ^ C[o,t](J^) '■ /(O) = x} with the subspace topology 
Tg^ = {U n Sx : U e r}. The set S* is a closed subset in Sx since the coordinate mapping 
/ !-;> f{T) is r-continuous. By the same property and the fact that the coordinate mapping is 
r-open, the r^^ -interior of S is int S = {f £ C[o^t]W ■ /(O) = x, f(T) > b}. We thus have 

inf Uf) = inf Ix(f) = inf inf Uf), and inf Uf) = inf inf 4(/). 

/ecl5 ^■" f€S ^■" a^bfGSa ^ ' feints ^■" a>b feSa ^ ' 

Using Proposition l2.2l and the fact that a^b > x{T), 



inf inf Ix{f) = inf inf 4(/) 



^7 ^ 1' 



a^bfeSa ^^-^ a>bfeSa ""'■" [l-e-27r]((j2/^) 

Consequently, S is an /a; -continuity set. D 

Now the decay rate under consideration can be determined. 

Proposition 2.4 Let b> x{T). Then 

[6 - (^ + (x - ^)e-T^)]2 
limelogF.(X^ ^ 5) = ^/_ ^,.,.j(>^^) ■ 

Moreover, the minimal cost path is as given in Proposition \2.2\ (with a replaced by b). 

Proof Apply 'Freidlin-Wentzell' to the I^; -continuity set 5: 

lim e log Fx (X^ ^ b) = lim e log P^ iX' G S) 

= -inf4(/) = -inf inf 4(/). 
/G5 a^bfeSa 

By the computations in the proof of Proposition l2.3i we obtain the desired result. The minimal 
cost path is directly obtained from Froposition l2.2[ D 

We mentioned in the introduction that there is an alternative method to compute the decay 
rate under study It follows relatively directly from the fact that X^ (with Xq = x) is normally 
distributed with mean ^t = x(T) = ^(1 — e~'^'^) + xe~"''^ and variance a'^{e) = ^(1 — e"^'^"^), 
in conjunction with the standard inequality ||TOl p. 19] 

1 1,,2 I 1+2 , 1 1,,2 



y + y ^ Jy y 

We have followed our sample-path approach, though, for two reasons: (i) the resulting most 
likely path is interesting in itself, as it gives insight into the behavior of the system conditional 



on the rare event, but, more importantly, (ii) it is useful when studying the counterpart of the 
decay rate for ROU (rather than OU), which we pursue in Section |3l 

We also note that 

(b--? 
lim lime log P^(X^ > b) = ^y— . 

It is known that the steady-state distribution of X^ with Xq = x is normally distributed with 
mean a/7 and variance e 0-^/(27). We conclude that this shows that the result is invariant 
under changing the orders of taking limits (T — )^ 00 and e — )■ 0). 



3 Transient asymptotics for reflected Ornstein-Uhlenbeck 

This section determines the decay rate ^ for ROU. For the moment we consider a setting more 
general than OU and ROU, namely stochastic differential equations with reflecting boundary 
conditions. Let D° G M be an open interval, dD and D denote its boundary and closure. Let 
u{x) denote the function giving the inward normal at x G dD, i.e. u{x) = 1 if x is a finite left 
endpoint of D and u{x) = —1 if x is a finite right endpoint of D. The reflected diffusion H"^ 
w.r.t. D is defined as the unique strong solution to 

dH^ = b{Hl)dt + yfeadBt + d^j , H^ = xeD, 

where \C% = Jq loDiHI)d\C^\s and CI = Jo ^iHs)d\C\s. Here \C% denotes the total variation of 
^■^ by time t. We assume that 6(-) is uniformly Lipschitz continuous and grows no faster than 
linearly, and o" is a nonzero constant. The existence and uniqueness of the strong solution is 
proved in |12|. 

Next we consider the solution to the deterministic Skorokhod problem for D: given a E 
C[o,oo)(l^)' there exists a unique pair (/i,/3) such that h G C'[o,oo)(-C'), and /3 G C[o,oo)(I^) of 
locally bounded variation that satisfy 



ht = at + Pu l/3|i=/ ldD{hs)dms, I3t= z^(/is)d|/3U 
Jo Jo 

We now recall the sample-path LDP for the reflected diffusion process, since it is consider- 
ably less known than the (standard) Freidlin-Wentzell theorem for the non-reflected case. We 
denote by H^ the nonnegative functions in H^ and by w a function from [0, T] to R. 

Theorem 3.1 (Doss and Priouret [3, Thm. 4.2]) Ifb{-) is uniformly Lipschitz continuous and bounded, 
and a is a nonzero constant, then H'^ satisfies the LDP in C[o,t](-D) '^^^h the rate function 

I{h) = inf -^2 / {K - b{ht) - v{ht)utl9D{ht)?'^t. 
ifhe H^ and 00 else. 



For reflected diffusions with a single reflecting boundary at 0, we can identify uj{t) and have 
the following explicit expression of the rate function. As usual, we define x+ = max{a:, 0} and 
x~ = — niin{x, 0}. 

Proposition 3.2 Let D = [0, oo). When 6(0) ^ 0, H^ satisfies the LDP in C[o,t]([0, oo)) with the 
rate function 



ifh € H^ and oo else. When 6(0) < 0, H'^ satisfies the LDP in C[o,t]([0, oo)) with the rate fiinction 
^+W = ^^ {K-b{ht)fdt-^h{Qfj^ l{0}(Mdi. 



ifhe H^ and oo else. In short, for h G H^ and 6(0) eRwe have 



I^ih) = ^f [K - Kht)f dt - i^iK^r? I Mo}{ht) dt. 
Proof In this case, D = [0, oo), dD = {0} and v{0) = 1. The rate function becomes 

I+(h)= mf -^ f ih'^-biht)-u;tl{o}iht)fdt. 

We minimize for each t separately under the integral. If h[ — b{ht) < 0, then wj = is optimal. 
If h[ — b{ht) ^ and ht > 0, then 1 rQ}(/it)aJi = 0, which means that any value of lo is optimal. 
If h[ — b{ht) ^ and ht = 0, then uit = h[ — b{ht) is optimal. Hence uj^ = {h[ — b{ht))^ is the 
optimizer. It gives the following explicit expression: 

1 t-T 



ii h £ H^ and oo else. For any h G Cp^T] ([0, oo)) which is differentiable a.e., note that h[ = Oii 
ht = 0. Then we have 

I^ih) = ^l {h't- b{ht) - ^o}{ht)b{Orf dt 

= ^l {h't-b{ht)fdt+^j^ ^^^^ht){b{Orfdt+^l^ l{0}(/it)6(0)-6(0)dt. 



/O ^" Jo " JQ 

When 6(0) ^ 0, the last two terms are zero, and for 6(0) < they sum to 

-^KO)'! l{o}{ht)dt. 
This completes our proof. D 
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Theorem 13.11 requires b{-) to be bounded, but this is a condition that ROU does not satisfy. A 
careful inspection of the proof in ||3l or fH, however, reveals that if a is a constant the bound- 
edness requirement for b{-) can be dropped. Specifically, for constant a, this proof can be 
modified in the sense that one can directly apply the contraction principle, which needs uni- 
form continuity of b{-) only. As a result, Proposition l3.2l is valid for ROU. Above we observed 
(i) that the most likely path for OU was non-negative (Proposition l2.2|) , (ii) the rate functions / 
and I~^ for OU and ROU are the same as long as their arguments are nonnegative paths on [0, T] 
(Proposition l3.2D . This suggests that the decay rates for OU and ROU (and the corresponding 
most likely paths) coincide. 

The idea is now that we find the decay rate (HJ for ROU by using the sample-path results that 
we derived in the previous section for OU. Recall that the zeroth-order approximation of OU is 
x{t) = a/7 + {x — 0/7)6^'''*. It is readily checked that x{t) > when the starting point x ^ 0. 
So we still assume b > x(T) in the decay rate for ROU. We define 5+ := {/ G C[o,r]([0, 00)) : 
/(O) = X, f{T) ^ 6}, corresponding to the rare event {y^ € S^}, so as to compute the decay 
rate ©; the set 5"+ is defined as {/ E C[o,t] ([0, 00)) : /(O) = x, f{T) = a}. Below we keep the 
notation F^ for probabilities of events in terms of V^ when this process starts in x. 



Theorem 3.3 Let b > x{T). Then, similar to the result of Proposition\ 

[b - {^ + {x - ^)e-^^r 



lime log P^(y^ ^ 6) 



.7 ■ \ 7' 



Moreover, the minimal cost path is as given in Proposition \2.2\ (with a replaced by b). 

Proof Since 6(0) = a > 0, by Proposition l3.2l V^ satisfies the sample path LDP in C[o,t]([0> o^)) 
w^ith the rate function 



It{h) = i^j {K-a + ^ht)fdt 



ii h £ H^ and 00 else. By an argument that is similar to the one used in the proof of Proposi- 
tion l2.3l 5+ is an /+ -continuity set. Then, 



lime log P^(y^ ^ 6) = limelogP^(y' g 5+) = - inf I+(h) = - inf inf /+(/i). 

e-)-0 e-)-0 ■ heS+ "■^bheS+ 



We have 

inf /+(/i) = mf\^ [ {h[ + ^ht-afdt,heH^nS+\ 

^ infl^ I {h[ + jht- afdt, he H^nSal 

The optimizer f* of inf f^Sa ^x{f) is always positive, for any starting point x ^ 0, due to Propo- 
sition |2!2l That is, /* G S+. Conclude that inf ^g^+ /+(/i) = inf /g^^ Ia;(/). Then the results 

10 



follows immediately from Proposition l2.2l and a ^ b > x{T), 

lim6logP.(y^^6) = -inf / .7 -2,ri72/ , 

[b - {^ + jx - ^)e-^^)? 
[l_e-27T](^2/^) • 

This proves the claim. D 

4 Transient asymptotics for doubly reflected Ornstein-Uhlenbeck 

This section computes the decay rate Q, but now for DROU. The case of DROU corresponds 
to choose the set D = [0, d] in Theorem 13.11 We can still derive an explicit expression for the 
optimal uj*, and hence have the following simplified rate function. 

Proposition 4.1 Given D = [0, d], the rate function I{h) in Theorem l3l] can be rewritten as follows: 

1 '■^ 



I^^ih) = ^J {K- b{ht) - l{o}(/it)6(0)- + l^,}{ht)b{d)+y dt 
ifh G H++ :={f eH^: 0< f <d}and 00 else. 
Proof Since dD = {0, d}, the rate function becomes 

I++{h) = mf ^ / {h[ - b{ht) - [l|o}(/it) - l{d}{ht)]ujtf dt. 

Let UJ* denote the optimizer of the above problem. We discuss the value of w* in two cases. 

• Case 1: h'-f. — b{ht) < 0. If ht G (0, d), then uj^ can be any value; if ht = 0, then uj* = 0; if 

ht = d, then uj* = — (/i^ — b{ht)). 

• Case 2: h[ — b{ht) ^ 0. If ht G (0, d), then uj^ can be any value; if ht = d, then cj^ = 0; if 
ht = 0, then uj* = h't — b{ht). 

As a consequence we have the following explicit expression: 

I^^{h) = ^ / {h't- b{ht) - l{,}{ht){h't - biht))+ + l{d}{ht){h't - b{ht))-fdt 

iih £ H++ and 00 else. Also, h G [0, d] and h is differentiable a.e. imply that Vi G (0, T), h't = 
when ht = or ht = d. So the above expression can be further simplified to 

I^^ih.) = ^l [K- Kht) - l{o}(/ii)("&(0))+ + l{d}{ht){~b{d))-f dt 
iih £ H^^ and 00 else. D 
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For DROU, 6(-) is bounded, and as a consequence it fulfills all requirements in Theorem 13.11 
and its rate function can be obtained by Proposition 14. 1 1 directly. Now, by a similar argument 
as employed in the last section, we prove that the decay rates (H) for OU and DROU (and the 
corresponding most likely paths) coincide. Recall from the introduction that we have assumed 
a/7 < d for DROU throughout this paper. Under this assumption, the zeroth-order approxi- 
mation x{t) belongs to (0, h) when the starting point 2; is in [0, d]. 
We consider crossing levels d ^ a ^ b > x{T). Define 

S++ := {/ e qo,T]{[0,d]) : /(O) = xJ{T) ^ 6}, 

so that our rare event corresponds to ZJ € 5++; the set S^"*" is defined as {/ G C[o,t]([0, d]) : 
/(O) = X, f{T) = a}. Finally, we arrive at the main result for DROU. 

Theorem 4.2 Letd^b> x{T). Then 

[b - {^ + {x - ^)e-^^)]' 



lime log P^.(Z|, ^b) 



1 ' 



Moreover, the minimal cost path is the one given in Proposition [ 



Proof Since 6(0) = a > and b{d) = a — ^d < 0,hy Proposition 14. 11 Z*^ satisfies the sample- 
path LDP in C[o,T]([Oi d]) with the rate function 

a h £ H^^ and 00 else. Due to Proposition I2.2[ it holds that /* = arginf/gg^ Ix{f) is in [0, d] 
on t G [0, T], when the starting point x G [0, d], a G [0, d] and a/7 < d, i.e., /* G S^~^. As an 
immediate consequence, /* = arginf ,gg++ /++((^), and inf^g^++ I^^{h) = infjgg^ Ix{f)- The 
rest of proof is similar to that of Theorem l3.3[ D 

As a side remark, we mention that if d < a/7, then the rate function reads 

1 l'^ 

It^{h) = ^ J {h't-a + lht + l{d}iht)ia - -fd))^ dt. 

In this case, the upper boundary does affect the transient behavior of DROU; this is in line with 
the intuition that in this case the process is 'pushed' towards the upper boundary. It is noted, 
however, that the functional w^ithin the integral is not differentiable with respect to ht, which 
seems to make analytically solving the associated variational problem challenging. 

5 Central limit theorem of loss process 

The main objective of this section is to derive a central limit theorem for the loss process Ut, 
for t large. We do so relying on martingale techniques. A similar procedure can be followed 
for the idleness process Lj. 
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Let /i be a twice continuously differentiable function on M, and Z be the DROU process defined 
earlier. By Ito's formula, we have: 

2 

Ah{Zt) = ((a - -iZt)h'{Zt) + y/i"(Zi))dt + ah'{Zt)dBt + h'{Zt)dLt - h\Zt)dUt. 
Based on the key properties of L and U, this reduces to 

dhiZt) = {Ch)iZt)dt + ah'{Zt)dBt + h'{0)dLt - h'{d)dUu (3) 

where the operator C is defined through 

dx 2 dx^ 
We first solve the following ODE with mixed conditions. 

Lemma 5.1 The ODE with real variable right hand side g G R 

{Ch){x) = q, ^ X ^ d, 
such that h{0) = 0, h'{0) = 0, and h'{d) = 1, has the unique solution 

a'_Wid) 2qu r rW{v) 

Q = Qu '■= 7, 1 h(x) = — TT / / T77r-:dv du, 

2 j;lw{v)dv ^ ' ^^ Jo Jo W{u) 

where 

W{v) := exp -^ j- 

Proof By applying reduction of order, the ODE can be written a system of first-order ordinary 
differential equations: 

h'f \ f/ \ f'f \ , 2a-27x 2q 

h (x) = / X , / (x) + 2 /^ = -^■ 

The integrating factor of the second first-order ODE is W{x). Hence, 

Ci_ 2q r W{u) 
W{x) ^a^ Jo W{x) 



fi^) = T7^ + -2l T777z(dn. 



Then the general solution is 

h{x) =C2+ r f{u)du =C2+Ci r TT^Au +% r r ^4^^ d«. 

Jo Jo W{u) (j2 70 70 W{u) 

Then the initial conditions /i(0) = 0, /i'(0) = uniquely determine the values of Ci, C2, while 
h'{d) = 1 uniquely determines qu- Hence, we obtained the desired unique solution. D 
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Proposition 5.2 The loss process U satisfies the central limit theorem, with rjjj defined in @, 

Ut - qut 



7V(0,??^), ast 



oo. 



Proof We insert the unique solution h{x) of Lemma ISH into ^. Since h'{0) = 0, h'{d) = 1, and 
{Ch){Zt) = qur we have the following integral expression: 



Ut - qut + h{Zt) - h{Zo) =a I h\Zs)dBs. 

Jo 



We then observe that Mt := Ut—qut+h{Zt)—h{Zo) is a zero-mean square integrable martingale. 
As usual, (M) denotes the quadratic variation process of M. By the ergodic theorem IH p. 134], 

t-i(M)t = t-V^ f h'iZsfds ^a^ f h'{x)\{dx) =: r/^, (4) 

Jo Jo 

where vr is the stationary distribution corresponding Zf. The density of vr is obtained in IIT4l 
Prop. 1] as, with N{m, s"^) denoting a Normal random variable with mean m and variance s^. 



d ^ / , / a cj^ , 



»«"'7'C)*^ 




f{('=-'^)S 



a\ / 27 



'^ (^-f)v^ -'^ (-f)v^ 



where c/? and <1> are the density function and cumulative density function of a standard Normal 
random variable. Then, 



Then by the martingale central limit theorem [8 . Thm. 2.1], t 2 Mt =^ AA(0, r/^) as t — )■ 00. Since 
Z G [0, d] and h is continuous, /i(^) is bounded. So, 

hiZt)-h{Zo) 

Vt 

a.s. as t — ;■ 00, which implies the claim. D 

The loss process at can be treated analogously. Define 

Proposition 5.3 The loss process L satisfies the central limit theorem 

:^L^M^AA(0,77i), flst^oo. 

Vt 
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6 Discussion on large deviations of loss process 

Realizing that qu is the mean rate at which U increases, it is seen that for c > qjj, the event 
{Ut > ct} for t large is rare. We analyze the decay rate of its large deviations probability 
limj_j.oo i^^ '^og¥z{Ut > ct) by appealing to the Gartner-Ellis theorem. Here and below we use 
the notation '^z{-) for the probabilities of events in terms of Z if Z starts in z, and "^zi") the 
corresponding expectation. 

To apply this result, our first task is to prove that t~^ logE^ ex.p{9Ut) converges, as t — )■ oo, to 
a function ip{9) (pointwise in G M; this function is often referred to as the limiting cumulant 
function). We apply Ito's formula to obtain a martingale which contains exp(6'C/t). Then we set 
h'{0) = 0, h'{d) = 9. Similarly to what we did above, we obtain 

9Ut - f {Ch){Zs)ds + h{Zt) - h{Zo) = a [ h'{Z,)dB,. 
Jo Jo 

The quadratic variation process of cr/p /i'(Zs)di?s is u^ Jq /i'(Zs)^ds. Since Z G [0,d] and h 
is twice continuously differentiable, Novikov's condition is satisfied. We have the following 
exponential martingale: 



Mt = exp (a I h'{Zs)dBs " y / h'iZ^fdi 



Also, 



/ ft (j2 

Mt = exp i9Ut - / [{Ch){Zs) + —h'{Zsf]ds + h{Zt) - h{Zo) 
Since h{Z) is bounded and E^Mj = 1 for all t ^ 0, 

^ log E, exp (9Ut - f ({Ch){Zs) + ^h'{Zs)A ds) ^0, as t ^ oo. (5) 

Suppose the integrand in (|5l) is merely a function of 9, i.e. {Ch){Zs) + ((7^/2) h'{Zs)'^ = 'ip{9), 
then we will have the desired limit, i.e., 

lim -logE^exp(6'C/t) = ^(6*). 

t— J-oo t 

We propose to prove this as follow^s. Firstly, by assuming the existence of ip{9), we are to solve 
the second-order nonlinear ODE 

2 

{Ch){x) + —h'{xf = ij{9), ^ X ^ d, (6) 

such that /i(0) = 0, h'{0) = 0. Then h'{d) is a function of ij{9), so we can write it as h'{d, ip{9)). 
Secondly, we try to prove that there exists a unique root ip{9) of equation h'{d, 'il){9)) = 9 for any 
G M. It is readily seen that Q with the given initial conditions has a unique solution. Letting 
h{x) be this solution, we can apply change of variables y{x) = exp(/i(x)), so as to obtain 

aV{x) + (-27a; + 2a)y'ix) - 2il,{9)y{x) = 0, y(0) = 1, y'(0) = 0. 
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We can solve it explicitly, but we cannot verify the existence of a function Tp{6) such that 
h'{d, ip{0)) = 9. This problem is caused by fact that the solution is rather involved (expressed 
in terms of Kummer's series [9]). 

In the context of reflected Brownian motion Zhang and Glynn fl^ managed to solve the prob- 
lem of identifying the limiting cumulant function ^(0) using the approach followed above. 
This is an example which shows that it is not alw^ays possible to apply methods that work in 
the case of constant drift in cases with state-dependent drift (such as DROU); cf. f\M . 
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